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We apply a recent quasiparticle model for the electronic properties of metallic helimagnets to 
calculate the transport properties of three-dimensional systems in the helically ordered phase. We 
focus on the ballistic regime t 2 Tcf >> 1 at weak disorder (large elastic mean-free time r) or 
intermediate temperature. In this regime, we find a leading temperature dependence of the electrical 
conductivity proportional to T. This is much stronger than either the Fermi-liquid contribution (T 2 ) 
or the contribution from helimagnon scattering in the clean limit (T 5 ^ 2 ). It is reminiscent of the 
behavior of non-magnetic two-dimensional metals, but the sign of the effect is opposite to that in 
the non-magnetic case. Experimental consequences of this result are discussed. 
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I. INTRODUCTION 

The electrical transport properties of metals have given 
rise to various surprises over the last thirty years. Within 
a nearly-free electron model with quenched or static dis- 
order, and to lowest order in the impurity concentration, 
the Boltzmann equation is exact and yields the familiar 
Drudc expression for the electrical conductivity, 

<t = ne 2 r/m e , (1.1) 

with n the electron number density, e the electron charge, 
m e the effective electron mass, and r the elastic mean- 
free time between collisions, which is weakly temperature 
dependent. Corrections to this result, in an expansion in 
the small parameter 1/epT, with e-p the Fermi energy, 
turned out to be very interesting. To analyze them, one 
needs to distinguish, in the thermodynamic limit, be- 
tween the diffusive regime at strong disorder or low tem- 
perature, Tt <C 1, and the ballistic regime at weak dis- 
order or intermediate temperature, Tt 1. (The latter 
regime should not be confused with ballistic transport in 
mesoscopic systems, where the mean-free path is large 
compared to the system size.) In three-dimensional (3D) 
simple metals in the diffusive regime, the leading correc- 
tion is nonanalytic in the temperature 

(Tt) 1 / 2 

Sa-wi, oc tJo —, \n~ (3D, diffusive) . (1.2a) 

Oft) 2 

The sign of this effect is positive, which reflects a neg- 
ative T-independent contribution that is non-universal 
(i.e., depends on an ultraviolet cutoff). The effect thus is 
localizing, i.e., it decreases the conductivity compared to 
the Drude value. In two-dimensional (2D) systems the 
effect is even more dramatic,— 

(5o- WL oc a ln ( rr ) (2D, diffusive) . (1.2b) 
e-pT 



These results have been reviewed in Ref. 0. The log- 
arithmic divergency in 2D perturbation theory signals 
a breakdown of transport theory, and the behavior at 
T = is insulating, albeit very weakly so. These ef- 
fects in non-interacting electron systems in general, and 
the logarithmic temperature dependence in 2D in par- 
ticular, are known as "weak localization" . They can be 
understood in terms of constructive interference in the 
electron-impurity scattering process,— or in terms of the 
exchange of certain soft or massless diffusive modes (ei- 
ther "Cooperons", or "diffusons") between electrons. 

Taking into account the screened Coulomb interaction 
between electrons leads to additional effects, and consid- 
erably enhances the complexity of the calculations. In the 
absence of quenched disorder, Fermi-liquid theory accu- 
rately describes the behavior and leads to a conductivity 
given by the Drude formula 

a c = ne 2 T C /m c , (1.3a) 

with a Coulomb scattering rate proportional to T 2 ,® 

TT 3 T 2 

l/rc = T — ■ (L3b) 

In the presence of both disorder and a Coulomb interac- 
tion, and in the diffusive regime, the leading corrections 
in the diffusive regime are qualitatively the same as those 
shown in Eqs. (|1.2[) . 7 ' 8 However, the physics behind the 
effects is different, and this is reflected in, e.g., the sensi- 
tivity of the results to an external magnetic field. These 
effects are often referred to as "Altshuler-Aronov effects" . 
In the ballistic regime, Tt 3> 1, the same problem has 
been analyzed for 2D systems by Zala et ali£ They found 

5a aa oc cr T/e F (2D, ballistic) . (1.4a) 

The sign of the correction depends on the interaction 
strength, but generically it is localizing, as it is in the 
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diffusive regime. More generally, the conductivity cor- 
rection can be written 

5a AA /a a = —f(T T ), (1.4b) 
e F r 

with f(x — > 0) cx In ir, and f(x — > oo) oc x in 2D. The 
crossover between the two limits was also determined in 
Ref.l and a physical interpretation in terms of scattering 
by Friedel oscillations was given. 

The most appropriate interpretation of <5oaa is as the 
result of a correction to the clean Fermi-liquid relaxation 
rate 1/rc- In the ballistic regime, this correction is small 
compared to 1/rc by a factor of 1/Tt <C 1. Adding this 
correction to 1/tc and 1/r according to Matthiesen's 
rule leads, in the regime where 1/rc *C 1/r, to Eq. 
(jl.4a[) . Another possible interpretation of <5oaa is as an 
interaction-induced temperature dependent correction to 
(To- In the Coulomb case, 5a aa is small compared to both 
ac and <jq, but we will see later that this is not neces- 
sarily true if the correction is mediated by a different 
interaction. 

A general explanation of the physics behind these non- 
analytic temperature dependencies is as follows. If there 
are massless excitations that couple to the relevant elec- 
tronic degrees of freedom, then the lack of a mass in the 
excitation propagator will lead to integrals that are singu- 
lar in the infrared, and these singularities are protected 
by a nonzero temperature T > (or frequency in the 
zero-temperature limit), which leads to a nonanalytic de- 
pendence on T. Notice that the massless excitation does 
not have to be of non-electronic origin: with a proper 
classification of modes they can themselves be electronic 
in nature without leading to double counting. For in- 
stance, in the case of the weak-localization singularities 
the relevant soft mode is either the Cooperon, or the dif- 
fuson, which are electronic particle-hole excitations that 
are distinct from, but couple to, the current mode whose 
correlations determine the conductivity. An example of 
a non-electronic mode that couples to the current is an 
acoustic phonon in the presence of an electron-phonon 
coupling. Since the latter is relatively weak, this leads 
to a T-dependence that is far subleading to the Fermi- 
liquid T 2 contribution, namely, the well-known Bloch- 
Griineisen T 5 law. 

In this paper we consider the exchange of a more exotic 
soft mode, namely, the hclimagnon excitation in metal- 
lic helimagnets. Helimagnets, the best known examples 
of which are MnSi and FeGe, form a class of magnetic 
materials with a preferred axis in spin space, charac- 
terized by a vector q. The magnetization displays fer- 
romagnetic order in any plane perpendicular to q, but 
the direction of the magnetization rotates as one moves 
along the q-axis, forming a spiral with wavelength 2n/q, 
where q = \q\ is the pitch wave number, see Fig. Q] A 
th eory of the ordered phase has been developed in Refs. 
[Toll 111 which we will refer to as paper I and paper II, 
respectively. The magnetic order splits the conduction 
band, as in the case of a ferromagnet, with a Stoner split- 




FIG. 1: Schematic magnetization pattern in a helimagnet. 

ting A that is proportional to the local magnetization. 
The Goldstone mode related to the broken symmetry in 
spin space, the helimagnon, has been studied in paper II. 
As one might expect, the helimagnon has an anisotropic 
frequency-momentum relation, 



wo(fc) = \jc z kl + c±_k\_ . (1.5a) 

The helimagnon is thus ferromagnon-like in the direction 
perpendicular to q, and antiferromagnon-like in the di- 
rection along q. Here and in what follows we take the 
pitch vector q to point in the z-direction, and k± is the 
component of the wave vector fc perpendicular to q. The 
elastic constants c z and cj_ are given by 

c z = J z \ 2 q 2 /4, (1.5b) 
c x = 7 ±A 2 /4. (1.5c) 

where hp is the Fermi wave number^ and j z and j± are 
numbers. For the model considered in papers I and II 
their values are 

72 = l/36 , 7x = l/96. (1.6) 

We will adopt these values for the purposes of this paper. 

For later reference we note that the helical magnetiza- 
tion structure is a result of the spin-orbit interaction, 
which is weak compared to the exchange interaction. 
Consequently, the pitch wave number q, which is pro- 
portional to the spin-orbit interaction, is small compared 
to the Fermi wave number, q/k-p <C 1. 

The metallic helimagnet MnSi in particular is a very 
well-studied material with many unusual properties that 
have been reported and discussed in detail in Refs. 
Il3lf3fl5l and papers I and II, among others. Here we 
focus on the electrical conductivity in the ordered phase, 
which displays helical order with a helix wavelength 
2ir/q w 180 A below a critical temperature T c ss 30 K 
at ambient pressure. Transport measurements in the or- 
dered phase have so far shown no significant deviations 
from Fermi-liquid T 2 -behavior. Consistent with this, a 
theoretical investigation of the clean limit in paper II 
showed that the leading effects of the helimagnons is a 
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T 5 / 2 correction to the Fermi-liquid behavior. Our moti- 
vation for investigating disorder corrections to the clean 
behavior is two- fold: First, the residual resistance of the 
cleanest MnSi samples puts them in the ballistic regime, 
and simple considerations suggest very interesting behav- 
ior in that regime. Second, MnSi shows very unusual 
transport behavior in the paramagnetic phase, namely, a 
T 3 / 2 -behavior of the resistivity over almost three decades 
in temperature^ The origin of this is not understood, 
but it is natural to speculate that remnants of the heli- 
cal order, which are observed in the same region, have 
something to do with it. It thus is prudent to first do a 
comprehensive study of effects of the helical order in the 
ordered phase, where conditions are more clearly defined. 

At a technical level, adding quenched disorder to the 
formalism used in papers I and II would be hard. We thus 
employ an effective model that was developed in Ref. [l6l . 
which we will refer to as paper III. Equations in papers 
I-III will be referenced in the format (I.x.y), etcJZ 

This paper is organized as follows. In Sec. [TT| we list 
our most important results for the convenience of readers 
who may not be interested in the technical details. In Sec. 
Mil we set up a transport theory based on the effective 
model, using the Kubo formalism. We first check, and 
demonstrate the efficiency of, the model by reproducing 
the clean-limit results of paper II in Sec. IIII Al and then 
proceed with the calculation in the ballistic limit in Sec. 
IIIIBI We discuss our results in Sec. [IVJ Some general 
points pertinent to transport theory are made in Ap- 
pendix [A] and several calculational details are relegated 
to additional appendices. 

II. RESULTS 

Since the details of the transport calculation for heli- 
magnets are quite technical, we first list our most perti- 
nent results without any derivation. 

In the helically ordered phase, the conductivity tensor 
is diagonal, but not isotropic. Taking the pitch vector 
q in z-direction, its nonzero elements are a zz = ctl, and 
<Jxx = Pyy = cr±- We find the leading correction to the 
conductivity in the ballistic regime in 3D to be propor- 
tional to tT. For a cubic crystal structure, as is the case 
for MnSi, we have 

— » ^(t) 2 (£)T <") 

where v is a parameter measuring deviations from a 
spherical Fermi surface, see Eq. (|3.2p below, and the pref- 
actor of the T-dependence is accurate to lowest order in 
the small parameter q/kp. 

This result is valid in a window of intermediate temper- 
atures. For asymptotically low temperatures, one finds 
diffusive, rather than ballistic, transport behavior, and 
for higher temperatures the ballistic behavior crosses over 
to either the Fermi-liquid behavior or the clean helimag- 
net conductivity, which is proportional to 1/T 5 / 2 , see 



paper II. For realistic parameter values (for known heli- 
magnets), the temperature window is 

T ba ii « T « Tbaa (e F r) (q/k F ) 3 (e F /A) 3 , (2.2) 

where Tbaii oc X/(epT) 2 is the lower limit of the ballistic 
regime. In this regime the dominant temperature depen- 
dence of the conductivity is given by Eq. (|2.ip . and the 
conductivity is 

(TL r L=tro + 8tr L ,x(T). (2.3) 

Notice that the sign of the correction is opposite to the 
Coulomb case, Eq. (|1.4a[) . That is, the effect of the heli- 
magnon exchange is antilocalizing. We will derive these 
results in Sec. IHII and discuss them in Sec. IIVI 

III. ELECTRICAL CONDUCTIVITY OF 
ITINERANT HELIMAGNETS 

We now set up a standard technical formalism for 
transport theory in the context of the effective model 
for metallic helimagnets that was given in Eq. (III. 2. 19). 
The electrical conductivity tensor <xy can be expressed in 
terms of an equilibrium current-current correlation func- 
tion by means of the Kubo formula 

i 

(JijliVL) = — - = 0)] , (3.1a) 

where 

Tr^ifi) = -e 2 T E 7/ E 3i(Pi)3i(P2) E 

xj)<T 2 (P2,«w 2 ))) ff2 (p2 I jw2 - ifi)). (3.1b) 

is the current-current susceptibility or polarization func- 
tion, with fj and 77 the fermionic fields. (. . .) denotes 
an average with respect to the action in Eq. (III. 2. 19). 
iuj = iuj n = i2iTT(n + 1/2) and ifi = i$l n — i2itTn 
(n = integer) are fermionic and bosonic Matsubara fre- 
quencies, respectively, and for simplicity we suppress the 
index n. j(p) = de p /dp is the current vertex, and for 
the electronic energy-momentum relation we use an ex- 
pression appropriate for a cubic crystal structure, 

with m e the electronic effective mass, and v — 0(1) a 
dimensionless coupling constant that measures deviations 
from a spherical Fermi surface. 

The conductivity as written above is actually the trans- 
port coefficient for the quasiparticles defined in paper 
III, which are described by the fermionic fields fj and rj. 
The physical conductivity is given in terms of the elec- 
tron fields tp and ip, which are related to the quasiparticle 
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FIG. 2: Graphic representation of the polarization function. 
The directed solid lines denote Green functions. 




equation 



G a (p,iuj) 



1 



(3.4) 



Here Go is the bare Green function, which is given by Eq. 
(III.2.10c). Equations ([331 E31> just shift the problem 
into the determination of the self energy E and the vertex 
function T. In order to evaluate the Kubo formula, it is 
most convenient to separately treat the cases with and 
without quenched disorder, respectively. 

A. Clean limit 




+ 




FIG. 3: Conserving approximation for the self energy and 
the vertex function. The dotted line denotes the effective 
potential. 



We first use the formalism developed so far to re-derive 
some of the results of paper II for the conductivity of 
helimagnets in the absence of any elastic scattering. This 
serves both as a check, and as a demonstration of how 
much simpler it is to evaluate the Kubo formula within 
the quasiparticle model compared to the model used in 
paper II. The calculation presented here is just a slight 
generalization of what is presented in Appendix ^ which 
serves to discuss the extent to which the approximations 
used are controlled. 



fields by the transformation given in Eqs. (III. 2. 9). How- 
ever, we will work to lowest order in the small parameter 
q/kp, and to this accuracy the quasiparticle conductivity 
is the same as the physical conductivity, as can readily 
be seen from Eqs. (III. 2. 9). 

The four-point fermionic correlation function in Eq. 
(|3.1b|) is conveniently expressed in terms of Green func- 
tions Q and a vector vertex function T with components 

r\ 

TTy(i^) = -e 2 T^2^^2ij l (p)g a (p,iuj)g a (p,iuj ~ id) 

xTl(p;iu,iuj-ifl), (3.3) 

see Fig. [2] This expression is valid if the Green function 
(Vcri (pij iui) Vo- 2 (P2, iw2)) is diagonal in both momentum 
and spin. For our effective model this is the case (whereas 
it was not the case in paper II), and Q is expressed in 
terms of the self energy E by means of the usual Dyson 



1. Conserving approximation for the conductivity 

It is well known that, in clean systems, care must be 
taken to treat the self energy E, which enters the Green 
function Q 1 and the vertex function T consistently in a 
conserving approximation ! 18 ' 19 The simplest consistent 
approximation, which is equivalent to the Boltzmann 
equation, is to treat the self energy in a self-consistent 
Born approximation, and the vertex function in a ladder 
approximation. These are graphically represented in Fig. 
[3] Analytically we have integral equations 

E^foiw) = -T^i^F(fc;p-fc,p;ift) 

iO k 

xQ a (p — k,iu> — id), (3.5a) 
for the self energy, and 



T a (p; iu),iw — id) = ij(p) — ^ -^^Z V(k;p — k,p; iil') Q a (p — k,iu; — id') Q a {p — k,iuj — iVl' — iQ) 

xT a (p — k: iuj — iSl', ilu — id' — id), (3.5b) 



V 



I 

for the vertex function. V is the effective potential from For completeness, we list the bare Green function from 
Eqs. (III.2.18). Eq. (III.2.10c), 

G ,a(p,iu) = -. , , (3.6a) 
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where 

(3.6b) 

with £ p = e p — ep, er = 1,2, and the effective potential 
from Eqs. (111.2.18)^ 

V(k;p : p';ifl) = V oX {k,ityl(k,p)i(-k,p'). (3.7a) 
Here 

V Q = X 2 (q 2 /8m 2 c ), (3.7b) 

and 

1 q 2 /3k 2 



2N F uj 2 {k) - (in) 2 



(3.7c) 



is the helimagnon susceptibility. For later reference, we 
also list its spectral function 

x"(k,u) — Imx(k, i£l — > u + iO) 

- — - — —JTX [&( u _ w o(fe)) - S{u + u>o(k))] 



12N F k F uj (k) 



(3.7d) 



Finally, 



7(fe,p) 



2A 



is a vertex function. This specifies all input parameters 
for the two coupled integral equations (|3.5[) . 

Since the two spin projections do not couple, we can 
restrict ourselves to one spin projection at a time, which 
effectively reduces the problem to one of spinless elec- 
trons. In what follows, we consider the contribution from 
the pole wi(p) and drop the spin label elsewhere. In the 
end, the contribution from the pole W2(p) can simply be 
added. 



2. Solution of the integral equations 

The coupled integral equations (|3.5p can now be solved 
by following a slight generalization of the procedure out- 
lined in Appendix |XJ The conductivity is still given by 
Eq. (|A5ap , and Eqs. (|A8[) remain valid. We thus have 
the single-particle relaxation rate given by 



(3.8a) 

with the zeroth moment of the potential spectrum given 
by 

- ^2 TH E ^i(jp)) V"(k-p; fc,p; «), 



T 



k,p 



(3.7e) with the potential V from Eqs. (|3.7p . This leads to 



(3.8b) 



r (e) 



-u 2 q 2 r°° 



8m 2 fc F 7V F J_ 



du 



u + e 
T 



:\p±-(k ± - P± )}x"(k-p,u). (3.9a) 



Here riB(s) = l/(e 2: — 1) and tif(^) = l/(e 2: + 1) are the Bose and Fermi distribution functions, respectively. Since 
the susceptibility x is s °ft at zero wave number, to leading order in the temperature this can be rewritten as 



r (e) 



u 2 q 2 



8m 2 fc F 7V F 



jTdu [n B (~)+n F (^)] ^5>M*))%i(p))*£ [*i-(p.L-*x)]V(fe-P,u). 



(3.9b) 

The same result is obtained from Eq. (III. 3. 5) by averaging l/r(p, e) over the Fermi surface. Evaluating the integral 
leads to 



r (e) 



1,024 x^/4 Mj£ J It J i^J 7 °( £ / 2T )' 



(3.9c) 



with T g = Ag 2 /6fc|. Here 







sin 2 93 cos 2 99 cos 2 (2^) 



(3.9d) 



with 



17, <p) = 1 + 2^[t7 2 + 2(1 - rf ) sin 2 9? cos 2 <p] + u 2 ^ + (1 + 2ry 2 - 3ry 4 ) sin 2 tp cos 2 <p] 



(3.9e) 
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We have normalized g such that g(y = 0) = 1. 70 is the 
n = member of a family of functions defined by 



ln{y) 



K, 



(2n+l)/2 



32 



drr ( 2n+ l)/2 ^ B(r) 



-n P (r + 2 2 /)+n F (r-2y)] ) (3.9f) 



with 



K„ = 4 jT^ sin" z = 2" +1 r 2 /r(/i + 1) , 

(3.9g) 

where T denotes the Gamma function. The same result is 
obtained by integrating Eq. (II.3.29d) or (III. 3. 6a) over 
the Fermi surface. For an explanation of the physical 
relevance of the temperature scale T q , see the discussion 
after Eq. fl~5) below. 

Similarly, Eqs. (|A8[) still hold, and we find a transport 
relaxation rate 



du 



Tip 



u + e 
T 



ri(e) = -JV P 
with 



(3.10a) 



k,p 



xV"(k-p;k,p;u). (3.10b) 



This agrees with Eq. (II. 3. 39b) after integration over the 
Fermi surface. Explicitly, we find 



512 x 6 5 / 4 V k F 



5/2 



7l (e/2T). 



(3.10c) 

Here 71 is given by Eq. (|3.9fj) with n = 1. The conduc- 
tivity is given by 



e 2 fc F AT F 



1 



6m 2 7-oo 4T cosh 2 (e/2T) 1*1 (e) ' 
which leads to a Drude formula 

a = ne 2 Ti/m c . 



(3.11a) 



(3.11b) 



Here n = k F /6ir 2 is the electron density per spin, since 
we are considering an effectively spinless problem, see the 
remark at the end of IIII A ll The transport relaxation 
rate is 



l/n 



iAg{v) 
C x 



q 



A r- ^ ^ 



A 



T 



Tn 



5/2 



(3.11c) 



in agreement with Eq. (II. 3. 40a). Here 
256 x 6 5 / 4 f°° dx 



5 7T 



cosh (ac)7i(x) 



186 . (3. lid) 



In deriving Eqs. (|3T0b|) - (f3TTTd|) we have assumed a 
spherical Fermi surface whenever doing so does not make 



the integral vanish. As a consequence, the ^-dependence 
of the prefactors in Eqs. (|3. 10cp and (|3.11d[) is exact, 
within our model, only to lowest order in v. 

We see that the current effective model reproduces the 
results of paper II for the clean case, and we have also 
calculated some of the prefactors that were not given 
explicitly in paper II. 



B. Ballistic limit 

We now add quenched disorder to the action, using 
the standard impurity model with an elastic relaxation 
time r described in paper III. We then need to distin- 
guish between the diffusive limit, where the relaxation 
rate 1/r is large compared to the frequency or temper- 
ature in appropriate units, and the ballistic one, where 
the opposite inequality holds. In the diffusive limit, it 
is well known that an infinite resummation of impurity 
diagrams is needed to work to a given order in the dis- 
order. In the ballistic limit, this is not the case, and a 
straightforward diagrammatic perturbative expansion in 
the number of impurity lines is possible. This yields im- 
purity corrections to the clean conductivity. For the case 
of electrons interacting via a screened Coulomb interac- 
tion, this has been investigated by Zala et al.^ and the 
development in the present case follows the same general 
lines. 

It is convenient to include the elastic relaxation rate in 
the bare Green function in a self-consistent Born approx- 
imation, see Fig. 0] That is, instead of the bare Green 
function Go, Eq. (|3.6a)l . we use 



G(p,iu) 



1 



iuj — ijJ\{p) + (i/2r)sgno; 



(3.12) 



Here we absorb the correction to the bare elastic relax- 
ation rate that was discussed in Sec. III. A of paper III 
in 1/r. In addition to using G instead of Go, diagrams 
must be decorated with explicit impurity lines, which di- 
agrammatically are denoted by dashed lines with crosses, 
and which carry a factor 

u Q = 1/2kN f t. (3.13) 

The Green function Q, Eq. (|3.4[) . can now be written 

1 



G(p,iu) 



G- 1 {p,iuj)-5Y<(p,iuj) ' 



(3.14) 



where the self energy 5T, does not contain the simple 
impurity self-energy that is incorporated in G. 



+ 



I * > 



FIG. 4: Defining equation for the Green function G (thick 
solid lines) in terms of the bare Green function Go (thin solid 
lines) and the impurity factor ito (dashed line with cross). 
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We are interested in the leading disorder correction to 
the clean resistivity calculated in paper II, and in the 
leading temperature dependence of that correction. To 
find this, it suffices to work to first order in both the 
disorder and the effective potential^ and we can expand 
the conductivity up to linear order in ST, and the vertex 
function T. From Eqs. (|3~Tal |3~5)) and (|3~T4")| we find 



the following expression for the static conductivity <xy = 
ReliniQ^o a.ij(ifl — ► fMO): 

<*H =<r$ ) +5<Tf j + 5af j , (3.15a) 

with 



r (0) 



Sal 



*4 



1 ^.h'P-Jj'lK 1 ' ^ 



V 



2T J 4tt cosh^(e/2T) 



G R (p,e) G A (p,e)-Re (G R (p,e)Y 



(3.15b) 



2T J 4tt cosh^(e/2T) 
1 f de 



2Rc 



(Gfl(p, e)) 2 G A (p, e) SE R (p, e) + (G R (p, e)) 3 ST R (p, e) 



(3.15c) 



1 



2T J 4tt cosh 2 (e/2T) 



Re 



G fl (p, e) Ga(p, e) r,(p; e + iO, e - iO) 



(G fl (p,e)f rj-Cpje + tCe + tO) 



(3.15d) 



To write Eqs. (|3. 15[) we have performed the Matsubara 
frequency sums and have introduced retarded and ad- 
vanced Green functions G r .a{p, e) = G(p, ioj — > e ± iO), 
and a retarded self energy ST, R (p, e) = S'S(p, iui — > e+iO). 
Diagrammatically, these contributions to the conductiv- 
ity are shown in Fig. [5] In evaluating these diagrams, 
we again make use of the small parameter q/kp -C 1. To 
lowest order in q/kp, in many cases the Green function 
G can be replaced by the free-electron Green function, 
which greatly simplifies the integrals. 

We further notice that the conductivity tensor is not 
isotropic, since the integrand depends on the helix pitch 
vector q. However, simple symmetry considerations show 
that it is still diagonal, with different components in the 
directions parallel and perpendicular to q, respectively, 

Saij = S tj [8 iz 5<7 L + (1 - S iz ) Sa±] . (3.16) 

The diagrams can be classified as follows. Diagram 



(0) in Fig. [5^a) represents . To lowest order in the 
disorder, and in q/ky, it yields the Drude conductivity, 
Eq. (TO}, 

4 ?) = % a ° t 1 + (Veft, (q/k F ) 2 )] • (3.17) 

Diagrams (i), (iii), (vii), and (ix) contribute to 5tr s , 
and the remaining diagrams contribute to <5<r r . Diagrams 

(1) and (ii) in Fig. [5^b) do not contain explicit impurity 
lines, and hence need to be evaluated to next-to-leading 
order in the disorder. The diagrams in Fig. [5Jc) contain 
an explicit impurity line, and evaluating them to leading 
order suffices. 

1. Diagrams without explicit impurity lines 

Let us first consider the diagrams (i) and (ii) . Standard 
techniques yield 



(a) 




(o) 





FIG. 5: Leading disorder corrections to the clean conductivity. Solid lines denote the Green function G, dotted lines denote 
the effective potential, and dashed lines with crosses denote the impurity factor wo- 



(i) _ -Vo 



5a) > = 



de 



4tt T J_ x cosh 2 (e/2T) V 



1 fu-e 



Re [J+ + -' + (k)-4+-'-(k)] 



(3.18a) 



6*W 



-Vo 1 



8tt 



de 



cosh 2 (e/2T) 
+n F 



Re [4-' + -(k) i+r + ~w 4~ ,++ ( k ) + 4~ ,++ ( k )] 



(3.18b) 



Here the J arc defined by convolutions of Green functions, 





-> + (k) = 


V 


k)G R (p) G R (p)G A (p) G R (p-k), 


(3.19a) 




-•-(fc) = 


ji{p) h{p)ii k ^P)i( k ^P - 

V 


k)G R (p) G R (p) G A {p) G A (p-k), 


(3.19b) 


4~ 




ji(P)jj(p)'V( k >P)'V( k >P- 


k) G R (p) G A (p) G R (p - k) G A (p - k), 


(3.19c) 


4~ 


■+-(*) = 


ji(p)jj( k )l( k ^P)^( k ^P- 


k) G R {p) G A (p) G R (p - k) G A (p - k), 


(3.19d) 



r 



where G RlA (p) — G R . A (p 7 e = 0). Other convolutions noting retarded and advanced Green functions, and the 
are defined analogously, with the upper ± indices dc- 



9 



comma separating them denoting the momentum struc- 
ture of the convolution. In writing Eqs. (|3.18[) we have 
neglected contributions from other convolutions of four 
Green functions that are easily shown to be of higher or- 
der in the disorder than the ones we kept. For instance, 
a complete expression for diagram (i) contains contribu- 
tions from J,J ++ '~ and J^ ++:+ , which are subleading 
in this sense. Also, a complete evaluation of the dia- 
grams yields nominal contributions proportional to x\ 
the Kramers-Kronig transform of x" ■ These vanish once 
the real part is taken, as is to be expected: by Fermi's 
golden rule, to first order in the interaction potential, 
the scattering cross-section and hence the conductivity 
depend only on the spectrum of the potential. Finally, 
we have used the fact that the internal frequencies u and 
e in Eqs. (|3.18[) scale as the temperature T. To find the 
leading temperature dependence, we therefore can drop 
the frequency dependence of the Green functions, and 
this is reflected in Eqs. (|3.19|) . 

To evaluate the integrals in Eqs. (|3.19[) we work to low- 
est order in q/kp. We further neglect A, since in our effec- 
tive single-spin-projection model it amounts (at q = 0) to 
just a shift of the Fermi energy. That is, we replace coi(p) 
in Eq. (|3.8b[) by £ p . We further use a nearly- free elec- 
tron expression for £ p , i.e., we put v = in Eq. (|3.2|) . By 
comparison with paper II, we see that this does not quali- 
tatively affect the result, see below. These simplifications 
lead in particular to ji(p) = pi/m c , and to lowest order in 
the disorder the integrals can be evaluated in the famil- 
iar approximation that replaces the integration over \p\ 
by a contour integration over £, p r^ which we will refer to 
as the AGD approximation. Power counting shows that 
the leading individual contributions to 8a are of order 
8a oc r 2 T 3 / 2 . This should be understood as the second 
term in an expansion of l/(l/r + T 3 / 2 ). (We recall that 
the clean single-particle relaxation rate is proportional to 
T 3 / 2 .) We know from paper II, and from Sec. lIII Al above. 
that these terms must cancel, and that the leading tem- 
perature dependence at 0(r 2 ) is r 2 T 5 / 2 . These terms, 
and higher ones in the diverging disorder expansion, must 
be resummed to yield l/(l/r + T 5 / 2 ), which is the inverse 
of the clean-limit transport rate added to the elastic scat- 
tering rate according to Matthiessen's rule. Corrections 
to these contributions are smaller by a factor of 1/rT 1 / 2 , 
which leads to conductivity corrections 8a oc tT. These 
all cancel by the same mechanism that leads to the can- 
cellation of the r 2 T 3 / 2 terms, and this can be seen with- 
out performing the integrals. Finally, the convolutions / 
in Eq. (|3. 18b|) are subleading in temperature compared 
to the J by power counting: at 0(t 2 ) they contribute 
to the clean-limit T 5 / 2 term, and the leading corrections 
are again small by a factor of 1/rT 1 / 2 , which leads to 
conductivity corrections 8a oc tT 2 . We thus obtain the 
following result: 

5af) + 8af = 0(t 2 T 5 / 2 ) + o(tT), (3.20) 

where o(x) denotes terms that are smaller than 0{x). 
The arguments leading to this conclusion are outlined in 



Appendix [B] The term of 0(t 2 ) was interpreted above, 
and we will not calculate the leading temperature depen- 
dence of the term of O(r) since we will find contributions 
of 0(tT) from other diagrams. 

2. Diagrams with explicit impurity lines 

We now turn to the diagrams in Fig. [5jc), which carry 
an explicit impurity line. Their contribution to the con- 
ductivity is of O(r), and it thus suffices to calculate them 
to leading order in the disorder. Before we do so, we 
identify the small parameter that controls our disorder 
expansion. As we point out in Appendix [B] the ex- 
pansion parameter for the convolutions J that appear 
in the integrand in Eqs. (|3 . 18[) is 8 = l/vpk±r, with 
vp = kp/m e the Fermi velocity. According to Eqs. (|1.5[) . 
the transverse wave number scales as the square root of 
the helimagnon frequency, which in turn scales as the 
temperature by virtue of Eqs. (|3.7c|) and (|3.18p . The 
small expansion parameter is thus 

(5=1/4x6 1/ V(£ft) 2 T/A, (3.21) 

and this will turn out to be true for the diagrams in 
Fig. (5](c) as well. This is different from the Coulomb 
case, where the small parameter that controls the ballistic 
regime is and it will be important for discussing 

the size of the ballistic regime in Sec. IIVI below. 

The diagrams in Fig. |5{c) all contain six Green func- 
tions that factorize into two sets of momentum convolu- 
tions containing n and 6— n Green functions, respectively, 
with n = 3 or n = 4. Diagrams (iii) - (vi) contain the 
(3, 3) partitions, whereas diagrams (vii) - (x) contain the 
(4, 2) partitions. The same power-counting arguments 
that we employed for diagrams (i) and (ii), and that are 
explained in Appendix [B] reveal the following: 

First, to lowest order in the small parameter (q/kp) (i.e., 
replacing the helimagnon Green functions by nearly-free 
electron Green functions), only the (3,3) partitions con- 
tribute to O(rT), whereas the (4,2) partitions are of 
higher order in the temperature. That is, 

Sa (vii)-(*) = o(tT) (3 22) 

to lowest order in (q/kp), and we will evaluate all other 
diagrams to lowest order in this small parameter as well. 
We will come back to what happens to higher order in 
(q/kp) in Sec. |mB3| below. 

Second, for the transverse conductivity correction 8a± 
only diagram (iii) contributes to 0(tT), whereas for 8a^ 
the other (3, 3) partitions also contribute. 

In addition, by considering the reality properties of the 
convolutions involved, one finds that, third, diagram (vi) 
is given in terms of the real part of a convolution that is 
purely imaginary, and hence does not contribute. There- 
fore, in order to obtain the transverse conductivity cor- 
rection 8a± to leading order in the small parameter 8, 
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Eq. (|3.21|) . one needs to calculate only diagram (iii). For 
the longitudinal correction Sa^ one needs to also consider 
diagrams (iv) and (v). 

Finally, a cursory inspection of the integrals in addition 
to power counting shows that the terms that contain 
a bosonic distribution function (in analogy to the first 
terms in Eqs. (|3.18a[) and (|3.18b[) . respectively), have a 
potential to be of 0(tT In 8) rather than of 0(tT). How- 
ever, the leading contribution to diagram (iii) does not 
contain such terms. Diagrams (iv) and (v) do, but the 
logarithmic terms cancel between these two diagrams, 
and this can be seen without performing the integrals. 
We thus conclude 



Sa± = 5a ( l H) +o{tT) = 0{tT)+o(tT), 

Sa L = Sag® + Sa ( ™' v) + o(rT) 
= 0(tT) + o(tT). 



(3.23a) 



(3.23b) 



a. Diagram (iii) After the above preliminary con- 
siderations, we now evaluate diagram (iii). The leading 
contribution can be written 



Sa 



{Hi) 



8?r T 



de 



1 



cosh 2 (e/2T) V 



k 



du 

IT 



XJlp 



T 



x"(k,u)Kj 



(fe). (3.24) 



Here 



Kt 



"(A;) = 



f^E^f G R { P )G R {p)G A {p) 



V 



v 



2ni k 2 N F 2 



+ 0(r), 



(3.25) 



and 



L++--(fc) 



V 



\ E 7(*; P) V - k) G R (p) G R (p) 
xG A {p - fe) 

(3.26) 



a 2tt N F m 2 c 2 
iv — .n.o +0(l/r, fc_|_) 



3 A 2 fc^ 



The second lines in Eqs. (|3.25[) and (|3.26p are easy to 
obtain in the AGD approximation. Only the term pro- 
portional to k± in 7(fe;p), Eq. (|3.7c[) . contributes to the 
leading temperature dependence, hence the proportion- 
ality to v 2 . We again have dropped the frequency depen- 
dence of the Green functions, since it does not contribute 
to the leading temperature dependence. Consequently, 
the integral over e in Eq. (|3.24[) can be performed. Us- 
ing the fact that the helimagnon spectrum x" is an odd 
function of the frequency, we can write 



Sa 



-u V 1 



4tt 



iS/ -x"(k,u)C(u/2T) 

(fe), (3.27a) 



xK++-L++-< 



with 



C(x) — cotha; — a:/sinh 2 x. 



(3.27b) 



We next cast the expressions corresponding to dia- 
grams (iv) and (v) in an analogous form, before perform- 
ing the final integrals. 

b. Diagrams (iv) and (v) Using the same techniques 
as for diagram (iii), we find for the leading contributions 
to diagrams (iv) and (v) 



Sa. 



(iv) + (v) 



-U V 1 



2vr V 

xM 



E 



du 



X "(k,u)C(u/2T) 



k •'- 00 
+- 



(fc)M+-'+(fc). 



(3.28) 



Here 

M+-'+(fe) 



^V-7ftp) G R (p) G A (p) G R (p - k) 
p 



3. T/ie conductivity in the ballistic limit 



(3.29) 



Before we collect our results, we return to the question 
of the diagrams that do not contribute to 0{tT) to low- 
est order in q/k-p. A calculation shows that at the next 
order in q/kp they do contribute, i.e., their behavior is 
the same as that of the diagrams we kept, only the pref- 
actor carries an additional factor of (q/k-p) 2 . For diagram 
(vii) this is demonstrated in Appendix [Cj for others, the 
results are analogous. We note that if one wanted to keep 
these terms, one would also have to take into account the 
difference between the quasiparticle conductivity and the 
physical conductivity that was mentioned at the begin- 
ning of Sec. Mil 

Collecting our results, we now have 



■in 



L++--{k) + 2M+-' + (fc) NC-' + (k 



(3.30) 



Here we show only terms that contribute to the leading 
temperature dependence of Sa. With help of Eqs. (|3.25l 
HOB] 15391 [337b) . and (|3~7d)) the final integrals are easily 
performed. We find 



Sa± 



576 



ao 





(t) 2 


' A 


T~ 


Gf) 




2e^ 





Sai, = 3Sa i 



(3.31a) 
(3.31b) 



This result is valid in a temperature regime Tb a n < T < 
T q , as explained below. A is an ultraviolet energy cutoff, 
which must be imposed, as in the Coulomb case^ since 
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only the hydrodynamic contributions to various parts of 
the integrands have been kept. The cutoff-dependent 
part of 8a is temperature independent; it is an interaction 
correction to the Drude conductivity. The temperature 
dependent part is independent of the cutoff. Notice that 
the constant contribution to da is positive, i.e., the effect 
of weak disorder in conjunction with helimagnons is an- 
tilocalizing. Accordingly, the temperature correction to 
the conductivity is negative. 

A necessary condition for this result to be valid is that 
the parameter 8, Eq. (|3.21|) . be small, 



T>T baU = A/16 V6(e F r) 2 



(3.32) 



The ballistic temperature scale Tbaii defined in this way 
marks the lower temperature limit of the ballistic regime. 
Another necessary condition is related to the fact the the 
helimagnon resonance frequency ujq has the form shown 
in Eq. (|1.5a[) only for wave numbers k < q. As was 
explained in paper II, this defines another temperature 
scale T q = Ag 2 /6fc|, and in order for Eqs. (|3.31|) to hold 
we must have T < T q . This also identifies the order 
of magnitude of the UV cutoff: A = 0(T q ). The tem- 
perature dependent contribution to 8a in Eqs. (|3.31|) is 
thus a small correction to the constant contribution. For 
k > q or, cquivalently, T > T q , the resonance frequency 
is u>g{k) oc ^fc\_ ft 2 , and hence the components of k scale 
as k z ~ k± ~ T 1 / 2 . Repeating the power counting argu- 
ments of Sec. IIIIB2I (see also Appendix [Bj , this yields, 
for the temperature-dependent part of 5a, 



8a±(T) oc <5<tl(T) cx — ao 



which is valid for T > T„ 



t) 



3/2 



1/2 



Eqs. l[3"3T|) times (T q /T) 



1/2 



(3.33) 

We note that this is just 
so in the regime T > T q , 



effectively a factor of \/T gets replaced by v/T^ 

We will discuss additional temperature scales, and the 
size of the ballistic regime, in Sec. IIVI below. 



IV. DISCUSSION AND CONCLUSION 

We now discuss our results. First, we give a detailed 
discussion of the range of validity of our results, and of 
the various temperature scales involved. We then give 
semi-quantitative estimates for the size of the predicted 
effects. In evaluating these estimates, one should keep in 
mind that the qualitative dependences on various param- 
eters are accurate, and so are the ratios of temperature 
scales etc., but the number- valued prefactors are model 
dependent and should not be taken too seriously. 



A. Temperature scales, and relaxation rates 

We start by giving an alternative derivation of the 
small parameter for the weak-disorder expansion, Eq. 



(I3.21[) . Consider the Green function G given in Eq. 
(|3.12[) . With k the soft wave number, and u the soft 
frequency, the fc-dependence of the various terms in the 
perturbation theory is given by its spectrum, 



G"{p-k,u) = -TrA(u-wiCp-fc)), 
with A a Lorentzian of the form 

1 1/2t 



A(x) 



7r x 2 + l/4r 2 



(4.1a) 



(4.1b) 



In the clean limit, r — * oo, A(a;) turns into a delta func- 
tion. To determine the relevant scale, we recall the scal- 
ing of frequency or temperature with the components of 
the wave vector. From Eq. (|1.5ap we have 



uj (k) 



cj_ ft j 



(4.2) 



With p on the Fermi surface, i.e., u>i(p) = 0, we have, 
from Eq. (|3.6bj) . wi(p-k) = -p± ■ fcx/m e + 0(k z ). We 
now scale u with T, p with ftp, and k\ with T/y/c±. 
Keeping only leading terms, we can write 



G"(p-k,u) 



-tS- 



8/2 



(p-k ± ) 2 + 8 2 /A 



(4.3) 



with p and k the scaled vectors p and k, respectively, 
and 8 from Eq. (|3.21[) . This confirms the role of 8 as the 
small parameter for the disorder expansion. 

8=1 defines the temperature scale that was denoted 
by Tbaii in Sec. IIII B~3l and that we list here again: 



Tball 



1 



A 



(4.4) 



16V6 (cft) 2 ' 

For a given disorder strength, this defines the lower tem- 
perature limit of the ballistic regime. 
A second relevant temperature scale is 

T q = Xq 2 /6k 2 , (4.5) 

which was introduced in paper II. As explained there, 
it is the energy scale related to the crossover from 
the anisotropic helimagnon spectrum to an isotropic 
ferromagnet-like spectrum with u>o(k) cx k 2 . T q > Tbaii 
provided e F T > (3/8^6) 1/2 k F /q « 0Ak F /q. With 
q/k? ~ 0.02, as is the case in MnSi, this means 6ft > 20, 
which always holds for good metals. The ballistic con- 
ductivity correction is then given by Eqs. (|3.31[) in the 
temperature window Tb a n < T < T q , and by Eq. (|3.33[) 
for T > T q . For T < Tb a n one has diffusive rather than 
ballistic transport behavior. 

At this point it is useful to cast our result for the con- 
ductivity correction in the form of a correction to the 
relaxation rate. The Drude conductivity plus the ballis- 
tic correction is a = ao + <5<7baii, which implies a total 
relaxation rate in the ballistic regime 



1 



Tball 



1 A _ foball 
T V °0 



1 TXV 1 



t 1,152 



A T„ 



(4.6) 
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Here we have absorbed the constant contribution to the 
ballistic rate in the Drude rate, and we have taken the 
correction to the longitudinal conductivity. To obtain the 
total transport rate we also need to add the clean-limit 
rate 1/ri from paper II or Sec. IIIlAl above, which in the 
current context is obtained from a ladder resummation 
of diagrams (i) and (ii) in Fig. [5] and taking the clean 
limit. For v = 1 we have g\(y = 1) « 0.2, and from Eq. 
fl3TTc| we find, 



1 

n 



50 



5/2 



(4.7) 



for a total transport rate in the ballistic regime 



1 

r tr 



1 

— 1 


1 


n 


T 


A 


(' 


50 


U'F 



1 



Tball 
8 



T 



5/2 



350 r 



e F T 



Here we have put v = 1 in Eq. (|4.6p and have approxi- 
mated the numerical prefactor. It also is illustrative to 
recall the clean-limit single-particle relaxation rate from 
paper II, which for generic wave vectors is (see also Eq. 



S.p. 

clean 



cx A 



3/2 



(4.9) 



Notice that the clean transport rate is smaller than the 
clean single-particle rate by a factor of T/A, as was shown 
in paper II, whereas the ballistic transport rate is quali- 
tatively the same as the ballistic single-particle rate, see 
Eqs. (14. 6|) and (III. 3. 12). That is, the cancellation mecha- 
nism between self-energy contributions and vertex correc- 
tions that is characteristic for clean transport problems 
(and also holds, e.g., in the electron-phonon scattering 
problem) is not operative in the presence of quenched dis- 
order. As a result, 1/Vban is small compared to V^ican 
by a factor of 5, Eq. (|3.2ip , as was to be expected, but it 
is not necessarily small compared to the clean transport 
rate l/r c i ean . Rather, the ballistic behavior will cross 
over to the clean behavior at a temperature 



T 



1-5/2 



0.05A/(e F T) 2/3 , 



(4.10) 



which provides a third relevant temperature scale. A 
fourth one is given by the temperature where the ballistic 
rate becomes equal to the clean Fermi- liquid rate 1 /tfl = 
T 2 /eF, which is 



Ti 



10" 



1-2 



(4.11) 



T q , Ti_5/ 2 , and T\_2 all provide upper limits for the 
regime where the conductivity correction is given by Eqs. 



(13.31[) . We thus conclude that the latter are valid in a 
temperature window given by 



< T < Min (T„ Ti_ B/ 2, T x _ 2 ) 



(4.12) 



Let us first compare 7]_g/2 with Ti_2 by writing the 
latter as 



Ti_2 « Ti_ 



1 



5/2 



50(e F r) 1 /3 \k F J \ A 



(4.13) 



If A = O(cf), then Ti_2 *C T\_ 5 /2- In a weak hclimagnet, 
where A/e F « q/kp, this is still true due to the small 
factor 1/50 (e F r) 1/3 < 1. 
Similarly, we can write 



200 e F T \k F J \ A 



(4.14) 



We again conclude that Ti_2 is the smaller of the two 
temperature scales provided that 



A > 1 



ep ~ 5(epT) 1 / 3 \k F 



<1 



1/3 



(4.15) 



As long as this condition is fulfilled, Ti_2 is the smallest 
of the three lower bounds. We now compare Ti_2 and 
Tbaii by writing 



Ti_2 thus is larger than Tbaii provided 



(4.16) 



(4.17) 



In order for the inequalities (|4.15|1 and (I4.17| to be com- 
patible, we must have 



6ft > &:f/<7, 



(4.18) 



which is also roughly the condition for T q > Tbaii- This 
is not a very stringent condition, and will generally be 
fulfilled in reasonably clean systems. 

We conclude that, if the inequalities (|4.18|) and (|4.17[) 
hold, the ballistic conductivity correction is given by Eqs. 
(|3.31[) in the temperature regime Tbaii < T < T\-%. For 
lower temperatures the behavior crosses over to diffusive 
transport, and for higher ones, to Fermi-liquid behavior. 
If (|4~Tg)) holds, but (|4"T7|) is violated, then the Fermi- 
liquid T 2 behavior will mask the ballistic T-dependence, 
and will have to be subtracted in order to observe the 
ballistic effect. 

We also need to remember that due to the broken ro- 
tational invariance in a solid-state system, there actually 
is a term proportional to k\ under the square root in Eq. 
(|1.5a[) . but it has a small prefactor. As was explained in 
papers I and II, this becomes relevant for temperatures 
below a scale T so — T q (q/kp) 2 . In this context we fur- 
ther need to come back to our discussion of the screening 
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of the effective interaction given in Eqs. (|3.7p , see Ref. 
I20I . As was shown in Eqs. (III. 2. 25), screening modifies 
the temperature scale T so to T so = T q (q/k F ) 2 (qv F /X) 2 . 
Requiring that this is temperature scale is smaller than 
Tbaii leads to one more constraint, namely, 

e F T < {k F /q)(\/qv F ). (4.19) 

Combined with (|4.18[) this leads to 

k F /q < e F T < (k F /q){X/qv F ). (4.20) 

as a necessary condition for the ballistic conductivity cor- 
rection to be given by Eqs. (|3.31[) . If the condition (|4.19[) 
is not fulfilled, then the lower temperature limit of the be- 
havior calculated above will be given by T so rather than 
by Tball. 

B. Quantitative predictions for experiments 

We now give some quantitative estimates, using pa- 
rameter values relevant for MnSi as follows (see paper 
I, and the references and discussion therein): k F — 
1.45 A" 1 , q/k F = 0.024, e F = 23, 000 K, m c = 4m , with 
m the free-electron mass, qv F ~ 1,000K. The value of 
A is uncertain; the large magnetic moment suggests that 
A is close to e F , but it is possible that A is smaller than 
e F by a factor of 40. This uncertainty in the value of 
the Stoner gap is a substantial impediment for making 
experimental predictions, especially since the theory de- 
pends quite strongly on whether qv F is larger or smaller 
than A. Our calculations are valid for qv F < A. 

The residual resistivity of the cleanest samples in Ref. 
[T5I was po w 0.33 /ifJcm, which corresponds to e F r sa 
1, 000. If A ~ e F , this is inside the disorder window given 
by Eq. (|4~20)) . and the condition (|4~T8l) is easily fulfilled. 
If A is substantially smaller than e F , then the second 
condition in Eq. (|4.20[) will be violated, and the lower 
limit of the ballistic regime as calculated above will be 
given by T so rather than by Tbaii- From Eqs. (14.41 14.5f> we 
see that Tba.11 ^ s smaller than T q by a factor of about 4,000, 
and from Eq. (|4. 1 5[) we see that Ti_2 is smaller than T q 
as long as A/e F > 0.005, or A > 150 K. Finally, from Eq. 
flUn]) it follows that Ti_ 2 > T ba ii as long as A/e F < 0.1. 
We conclude that the parameter values of MnSi provide a 
sizeable ballistic regime. However, depending on the size 
of A, it may be necessary to subtract the Fermi-liquid T 2 
contribution to the conductivity in order to observe the 
ballistic correction. 

The absolute size of the effect, on the other hand, 
is very small. From Eq. (|4.1ip we estimate that Ti_2 
is at best, for the smallest conceivable value of A (1=3 
500K), in the mK range, and from Eqs. (|3.31[) we have 
\da h /(7o\ w 2 x 10- 7 (e F /A) 2 T/e F . For A = 500K) this 
yields I&tl/ooI ~ 4 x 10~ 4 T/e F . For temperatures on 
the order of T1-2, this makes for an extremely small ef- 
fect, and even at T 1 K, which requires subtraction of 
the Fermi-liquid contribution, the effect is small. 



It is conceivable that in other materials the effect is 
larger, or that artificial systems can be constructed, e.g., 
optical lattices, that have parameter values leading to a 
larger effect. The most efficient way to increase the ef- 
fect would be a larger helix pitch wave number. In real 
systems, the basic reason for the small absolute value of 
the effect is the prefactor (q/k F ) 3 in Eqs. (|3.3ip . This 
in turn reflects the fact that any effect of the helix will 
reflect the large size (on an atomic scale) of the helix, 
which leads to correspondingly small energy scales. The 
same comment holds for the conductivity in the clean 
case, see Eqs. (|3.11bl I3.11c[) . By contrast, if one writes 
the observed resistivity of MnSi in the disordered phase 
as p = Qq\1 + const, x (T/e F ) 3 / 2 ], then the experiment in 
Ref. [la yields const. = O(10 6 ). The anomalous tempera- 
ture dependence of the resistivity in MnSi is thus a very 
large effect that must be related to effects on small length 
scales, and is not likely to be associated with remnants 
of helical order. 



C. Conclusion 

In conclusion, we have applied the effective model for 
helimagnets that was derived in paper III to determine 
the effects of helical magnetic order on the electrical con- 
ductivity. In the clean limit, we reproduce the results 
obtained earlier in paper II, but the effective model al- 
lows for a much simpler calculation. We have applied 
this theory to determine the conductivity in the ballis- 
tic regime, which in helimagnets is characterized by the 
requirement y/ (e F r) 2 T/A 3> 1. Remarkably, we have 
found that the temperature correction to the resistivity 
in bulk helimagnets is linear in T, as it is in 2-d nonmag- 
netic metals. This analogy between 3-d and 2-d systems 
is a consequence of the anisotropic dispersion relation of 
the helical Goldstone mode or helimagnon. The abso- 
lute value of the effect, with parameters values appro- 
priate for known helimagnets, is very small due to the 
large size of the helix. We finally mention that the trans- 
port properties of helimagnets in the diffusive regime, 
\J (e F r) 2 T/A <C 1, remain to be investigated. Prelim- 
inary results suggest that they are less exotic than in 
the ballistic regime, with no effective reduction of the 
dimensionality. 23 
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APPENDIX A: ELECTRICAL RESISTIVITY DUE 
TO A GENERIC POTENTIAL 

For pedagogical reasons, and to make a few technical 
points that are not emphasized in the elementary litera- 
ture, let us consider the resistivity of nonmagnetic, spin- 
less electrons due to scattering by an effective dynamical 
potential. The familiar example of electron-phonon scat- 
tering, which leads to the Bloch-Griineisen T 5 law, is a 
particular realization of this generic case. The develop- 
ment in Sec. lIII Al follows the same logic. The only differ- 
ences are that the resonance frequency of the bare Green 
function is different, and that the potential depends on 
all three momenta in the scattering process, not just on 
the net transferred momentum. 

Consider spinless electrons interacting with a spin- 
independent, frequency-dependent, effective potential 
V(k,id). We assume that the spectrum of the poten- 



tial, V"(k,u) = ImV(fc,iQ ->■ u + iO), is soft at fc = 
and u = 0. The Green function is diagonal in both spin 
and momentum, 

S(p, iu) = -. c * . r , (Al) 

IU) — £ p — 2j(p, VjJ) 

where £ p = e p — fi, with e p the electronic energy- 
momentum relation, and /i the chemical potential. The 
self-consistent Born equation for the self energy E, de- 
picted graphically in Fig. [3l reads 

E(p, tw) = -^J2 T J2 V ^ in ) g (P id), 

p ifl 

and the integral equation for the vertex function T in a 
ladder approximation, also shown in Fig. [3j is 



r(p; iw, iu - id) = ij{p)~ — ^T^V{k,id')g{p - k,iu - id')g{p - k,iw - id' - id) 

k i£~t' 

xT(jp -k;iu -id' ,iw — id' -id). (A3a) 

If we define a scalar vertex function 7 by T(p; iuj, iuj — id) = ij(p) 7(p; iu, iu — id), we find that 7 obeys an integral 
equation 

y(p;iw,iu-id) = 1 - — V T V ' £f P ~ ^ V(k, id') Q(v-k,iw- id') Q(p-k,iu- id' - id) 

J (p) 

xj(p — fc; iuj — id', iu — id' — id). (A3b) 
I 



The polarization function and conductivity tensors are 
diagonal, atj(id) — 5ij<j(id). The sum over fermionic 
Matsubara frequencies in Eq. (|A3|) can be transformed 
into an integral along the real axis by standard methods. 
This procedure yields two terms where the frequency ar- 
guments of the Green functions lie on the same side of 
the real axis, and two other terms where they lie on op- 
posite sides. Only the latter contribute to the leading 
result as the self energy goes to zero. Since the real 
part of the self energy just renormalizes the Fermi en- 
ergy, and the imaginary part, which gives the relaxation 
rate, indeed goes to zero as T — > 0, we need to keep only 
these retarded-advanced combinations for the purpose of 
determining the leading low-temperature dependence of 
the conductivity. The Kubo formula for the static con- 
ductivity a = HmQ_>o Rcer(ii7 — > d + iO) becomes 

a = f. r de 1 Ly f -( i) 2 

3^ AT cosh 2 (e/2T) V ^ K3KP) > 

x \Q{p, e + i0)\ 2 j(jp; e + i0,e- iO). (A4) 
The Green functions in Eq. (|A4[) ensure that the dom- 



inant contribution to the sum over wave vectors in the 
limit of a vanishing self energy comes from p such that 
£ p = e. Furthermore, since e scales with T, for the 
leading temperature dependence we can neglect all e- 
dependencies that do not occur in the form e/T. In a 
nearly-free electron model, with a spherical Fermi sur- 
face with Fermi wave number kp, and j(p) = p/m Q with 
m e the effective electron mass, we thus have 

e»*g /~ & 1 A(e) 

3m 2 J_ x AT cosh 2 (e/2T) r (e) ' V ; 

Here we have defined 

A W = ^ 7(p;e + *0,e-i0), (A5b) 

p 

r (e) = ^$^(£ p )ImE(p,e + iO), (A5c) 
p 

and we neglect the real part of the self energy, which only 
redefines the zero of energy. 

Using analogous arguments, we find from Eq. (|A3b[) 
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that A obeys an integral equation 
A(e) = 1-N F [ duV"{u) 



T 



n-p 



T 



xA(e + M)/r (e + w). 



(A6) 



Here iib(i) = l/(e x — 1) and uf(x) = l/(e x + 1) are the 
Bose and Fermi distribution function, respectively, and 



V"(u) = 



1 1 

Si V 2 



Y,S(Ck)S^ p )V"(k-p,u)k-p/k 2 



fe,p 



(A7a) 

with Sf — 0-/V) Y^,k ^(£fe)> is an ^ = 1 average of the 
spectrum of the potential over the Fermi surface. For the 
purpose of finding the leading temperature dependence 
of the conductivity, it can be written 



V"{u) = V£{u)-V{'{u), 



(A7b) 



with 



2/V,.- 



= w I d PP (p 2 /2k F ) n V"(p,u). (A7c) 



The integral equation (|A6j) is not easy to solve. How- 
ever, in an approximation that replaces A(e+u)/To(e+u) 
under the integral by A(e)/ro(e), it turns into an alge- 
braic equation whose solution is 



A(e)=r (e)/r 1 (e), 



(A8a) 



Here we have used the fact that, in the limit of a small 
self energy, one find from Eqs. (|A5cP and (|A2[) . 



r (e) = -AT F / duV "{u) 

and we have defined 
ri(e) = -N F I duV"(u) 



(?) 



fU 



TIF 



TIF 



T 



(A8b) 



T 



(A8c) 

We see that the vertex function A effective replaces 
the single-particle relaxation rate T with the transport 
relaxation rate IV To see the relation between the two, 
we recall that the frequency u scales with the temper- 
ature. For potentials where the frequency scales with 
some (positive) power of the wave number, Ti will thus 



J 



depend on a higher power of the temperature as T — > 
than r. As an example, consider the case of elec- 
tron scattering by acoustic phonons, where V" (p, u) cx 
cp [S(u— cp) — S(u+cp)], with c the speed of sound. In this 
case, r (e) cx T 3 j {e/T), whereas Ti(e) cx T 5 7 i(e/T), 
where 



7n(y) = / dxx 2{n+1) [2n B {x) + n F (x + y)+n F (x- y)] 
Jo 

(A9) 

In this case, the single-particle scattering rate shows a 
T 3 dependence, whereas the transport scattering rate, 
and hence the resistivity, display the familiar Bloch- 
Griineisen law, a cx T 5 . 

Most of the technical development sketched above can 
be found in textbooks.— What is usually not stressed 
is the fact that the approximate solution, Eq. (|A8a[) , of 
the integral equation, (|A6p , yields the asymptotically ex- 
act temperature dependence (although not the prefactor) 
of the conductivity. The fact that it does has, to our 
knowledge, never been established within diagrammatic 
many-body theory (and it is not proven by the above ar- 
guments), but it can been seen from the fact that the 
asymptotic solution reproduces the lowest-order varia- 
tional solution of the Boltzmann equation^ The relation 
between a diagrammatic evaluation of the Kubo formula 
and solutions of the Boltzmann equation is complex, and 
will be discussed in more detail elsewherei^ 3 . 



APPENDIX B: POWER COUNTING FOR 
DIAGRAMS (i) AND (ii) 



Here we provide the arguments that lead to Eq. (|3.20[) . 
We first do a power-counting analysis of Eqs. (|3.18l) . 
From Eqs. lfOa| l3~7cl l3~18l) we see that the soft hc- 
limagnon wave number k scales with temperature as 
k z ~ kj_ ~ T. The frequencies scale as u ~ e ~ T, and 
x"(k,u) ~ 1/T 2 . Consequently, the conductivity correc- 
tions Sa^'**) scale as 5a ~ TJ for a given integrand J{k) 
(or 1(h)). 

First consider the integral J^ + ~' + {k), Eq. (|3.19a|) . 
For power-counting purposes, the integration variable p 
scales as T°, and the leading term in the vertex 7 scales 
as 7(fc,p) ~ k±_ ~ T 1 / 2 . A representation that suffices 
for power counting is thus 



J, 



++-,+ 



(k) cx fcj 



d£ / dr) 



1 



1 



1 



1 1 (^-i/2r) 2 Z + i/2r £-i/2T~v F kri 
I 



r 2 k 2 f V¥kT dx 



1 + x 2 



(Bla) 



in the AGD approximation. For VFk ^> 1/t we thus have factor of l/vphr ~ 1/rT 1 / 2 , or 



J, 



++-■+ 



(fc) cx t kj_/k, with corrections carrying an extra 



J2 



- + {k)~T 2 T 1 ' 2 +T 1 



(Bib) 
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which leads to da oc 

T 2 T 3/2 + tT ^ Analo 

gous arguments 

yield 



4 + - 


<~{k) ~ T 2 T^ 2 +T, 


(B2) 




f -(fc) ~ T 2 T l ' 2 +T, 


(B3) 




f+ (fc) ~ r. 


(B4) 



The convolutions /, compared to the corresponding </, 
carry an additional factor of k± ~ T 1 / 2 . In addition, 



the resulting vector nature of the integrand leads to an 
another factor of either kj_ ~ T 1 / 2 , or A: 2 ~ T. There- 
fore, the I carry an additional factor of T compared to 
the corresponding J. Terms that were dropped in writ- 
ing Eqs. (l3~18l) involved J+++<+, J+++ ~, J++-++, and 
, which are of higher order in the disorder by at 
least three powers of 1/r. Including terms of 0(tT), we 
thus can write the conductivity correction, Eqs. (|3.18|) . 



-Vf, 1 



de 



4tt rj. M cosh 2 (e/T) 
1 



«B (~ ) Re 



7-++-,+ 



( fc ) + 2 ^ 



+ -,+- 



(fc) 



■ Tip 



Re [jJ+-' + (*0 + 4"' + -(fe) - 4 + —{k) J+- ++ (fc)] 



2 V T 

The J can be simplified by means of partial fraction decompositions. For the relevant combinations one finds 



Re 



(k) 



2 l J 



(fc) 



cxr 2 fc1 ~ T 2 T 3/2 , 



Re [J+ + -' + (fc) + - 4 + -'-{k) - J+-' ++ (fc)] = o(tT°). 

I 



(B5) 

(B6a) 
(B6b) 



This leads to Eq. (|3~20|) . 



APPENDIX C: DIAGRAM (vii) 

Here we consider diagram (vii) in Fig. [UJc) as a pro- 
totype of a class of diagrams that do not contribute to 



the leading behavior of the conductivity if evaluated to 
lowest order in q/kp. The leading contribution to the 
conductivity correction from this diagram can be writ- 
ten 



J 



6a 



-u V 1 



de 



irm 2 T J^oo cosh 2 (e/2T) J-^ tt 



du ( u — e 



lm lpJ2 XR ( k > u } "^5Z 7 ( fe ' p ) G k(p) G A {p-k) 



7 7 ( fc ' p ') Gi? ( p ') G «( p ') Ga{ p'^ Ga{ p' - fe )' 



p' 



(Cl) 



which shows the (2,4) structure mentioned in Sec. MI B~2l 
The bosonic distribution function does not contribute to 
this diagram, so it can be at most of 0(tT). With the 
convolutions evaluated for q — 0, power counting shows 
that it is of 0(tT 2 ), and an explicit calculation confirms 
this. Now we expand the resonance frequency u>x(p), Eq. 
(13.6b[) , to first order in q: w\(jp) = £ p +p- q/2m c + 0(q 2 ). 
For the leading contribution to the first convolution in 



Eq. (jCTj) we then find 

7V+<-(fc) = -^5>(fc,P)G«(p)G A (p-fc) 



Xk 2 



k 1 k z 



+ 



kjq 
k 2 k F 



0(q 2 ) 



(C2) 



We see that, at linear order in q, a factor that used to be 
kjk ~ k z /k ± ~ T 1 ' 2 gets replaced by q/k F ~ T°. The 
same holds for the other convolution. As a result, the 
diagram is of O(rT), and an explicit calculation shows 
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that the dependence of the prefactor on ep/A and q/k-p 
are the same as for diagram (hi), with the exception of 
the additional factor of (q/k-p) 2 - We thus have 
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